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The electromagnetic field solution for a spherical invisibility cloak with an active device inside is 
established. Extraordinary electric and magnetic surface voltages are induced at the inner boundary 
... of a spherical cloak, which prevent electromagnetic waves from going out. The phase and handness 

00 ' of polarized waves obliquely incident on such boundaries is kept in the reflected waves. The surface 

, voltages due to an electric dipole inside the concealed region are found equal to the auxiliary scalar 

■ potentials at the inner boundary, which consequently gain physical counterparts in this case. 

CN ■ 

O ! PACS numbers: 41.20.Jb, 42. 25. Fx 

D ■ 

Pendry et al. pj firstly proposed the invisibihty cloak based on coordinate transforms where a "hole" is created in 
the transformed coordinate system and an object in the "hole" can be concealed from detection. Another method was 
' ' reported to produce similar effects in the geometric limit 0- Designs of invisibility cloaks based on other methods 
— ,' have also been published From the viewpoint of transformation method pj], the hole creation docs not 

' result in an electromagnetic vacuum but rather a complete separation of electromagnetic domains into a cloaked 
' region and those outside [11,01 • More precisely, a "true" cloak should not only cloak passive objects from incoming 
, waves, but also cloak active devices by preventing waves from going out and being detected. The effectiveness of a 
Q ' transformation based cloak design for hiding a passive object was first confirmed computationally in the geometric 
^ , optics limit in full- wave finite-element simulations (olliot. and in full- wave analytic scattering models with deeper 

^ ' physical interpretations pTi . [T^ . It has also been demonstrated experimentally with a simplified practical model . 
' ^ \ On the other hand, however, the electromagnetic wave behavior in this concealed region with an active device inside 
remains unknown, due to the fact that the concealed region, or the "hole" , created by the transformation method 
f~| ' does not exist before transformation and has no counterpart in the original Euclidian space. A rigorous mathematical 
Qh< treatment of cloaking shows that "finite energy solutions" to Maxwell's equations do not exist in the presence of active 
' ' sources inside the concealed region [l^ , which strongly implies that some special physical phenomenon must exist but 
. has not been revealed. Also, the exact electromagnetic field solution to the Maxwell equations in this case, as well as 
J> ■ whether waves can go out of the concealed region, still haven't been established or confirmed. 

T-H ] In this Letter, the exact field solution to the Maxwell equations with an active source inside the concealed region 
i: of a spherical cloak is established. We show that, electric and magnetic surface voltages are induced due to an 
] infinite polarization of the material at the inner boundary, which prevent the electromagnetic waves from going out. 
. The special property of the material at the inner boundary of the cloak, is able to keep not only the handness 

■ of polarization, but also the phase information, in the waves reflected from the inner boundary, whose behavior is 
different from a perfect electric/magnetic conductor (PEC/PMC), or the artificial soft and hard surface (SHS) in 
electrical engineering p^ . The induced surface electric and magnetic voltages are shown to be exactly equal to 
the auxiliary scalar electric and magnetic potentials at the inner boundary, respectively. Therefore, these auxiliary 
potentials, which are commonly introduced in accompany with vector potentials as strictly mathematical tools for 
most engineers [l6j . have physical counterparts at the inner boundary in this case. 

Figure[l]shows the configuration of a spherical cloak with outer radius i?2 and inner radius Ri . The cloak layer within 
Ci ' i?i < r < i?2 is a specified anisotropic and inhomogeneous medium with permittivity tensor e — erT^r + etOO + et(t>4> and 
permeability tensor /2 = iir'f"r+^t00+fit4'4'- According to Ref. it is chosen such that et/eo = /^t/Mo = R2/ {R2 — Ri) 
and e^/et = li-r/ li-t = {r — Ri)^ /r^ . Without loss of generality, we assume that the background material in the region 
r < Ri has permittivity ei and permeability ^1. A time-harmonic electric dipole is put inside as an active device. 
The electromagnetic waves from the dipole as well as the response of the surrounding environment can be decomposed 
into TE and TM modes with respect to f, corresponding to scalar potentials ^'te and ^tm, whose expressions in 
the current case have been derived in Ref. [ll|. 

Since TE and TM modes in such a radially inhomogeneous medium can be shown to be decoupled the 
derivations for these two kinds of modes are identical to each other. We start with the case that an outgoing TM 



'Electronic address; lchenhs@ewt.mit.edul 



FIG. 1: (Color online) Reflection and transmission of electromagnetic waves with an electric dipole inside the concealed region 
of a spherical cloak. Region 0: concealed region; Region 1: cloak layer; Region 2: free space. 



wave is excited in the concealed region. This outgoing wave will induce a standing wave in re gion r < i?i , an outgoing 
wave and a standing wave in region Ri < r < i?2, and an outgoing wave in region r > R2 [l7|, as shown in Fig. [1] 
Thus the scalar potentials in these three regions can be written as 

nil = lUkir) + R^''Mkir)]P,Ticose)e^"''f> (1) 

*TM = [d^Uniktir ~ Ri)) + f^'xnihir ~ Rl))] 

■ P^\cos0)e""^ (2) 

5-°"*, = r™Cn(fcor)Pr(cos0)e™^ (3) 

where ipn, Xn, and C„ are Riccati-Bessel Functions of the first, the second, and the third kind respectively; R^^^ ^ d*^, 
/*^, and T"^*^ arc the unknown expansion coefficients. Especially, R^^^ and T-^^^ arc called the general reflection 
coefficient and general transmission coefficient, respectively 

For the sake of illustration, the inner boundary of the cloak is set at r = i?i + (5 instead of r = i?i and then the 
limit 5 ^ is taken [l^ . Consequently, four boundary equations can be listed utilizing the continuities of tangential 
E and tangential H at the outer boundary (Eq.(|4]) and ([5])) and at the inner boundary (Eq.® and ([7])): 

i/^t[dt'i^'MR2 - Ri)) + f^^XniHR2 - Ri))] 

= l/^JI^T^''C{koR2) (4) 

l/fitKMkt{R2 - Rl)) + Xn(fct(i?2 - Rl))] 

= l/fioT^''UkoR2) (5) 
l/VJ^[Ciki{Ri + 5)) + + S))] 

= l/^t[d^^UktS) + fJi'x'niktS)] (6) 
l/Aii[C„(fci(i?i + 6)) + RT''i'n{ki{Ri + 6))] 

= i/Md'JMhS) + f^'xnihS)] (7) 

After solving all the equations, it can be obtained that d^^ — f^'^ = T'^^' = 0, indicating that no field exists in the 
cloak layer as well as the outside space. Meanwhile we can get that R^^^ ~ ~'(n{kiRi)/^pn{kiRi), which is important 
for later use. In addition, it follows that in the limit (5 — > 0, fn^ Xn{kt5) is nonzero. Obviously, the value of this 
product in Eq.® is nonzero only at the inner boundary, bringing out the discontinuity of the tangential E field across 
the inner boundary. This result is very interesting. Since both /.i and e are finite everywhere and no conductive media 
exist, it is unreasonable to include surface current to support this discontinuity, as the common boundary conditions 
do. It should be noted that due to the same reason the displacement surface currents introduced in the cylindrical 
cloak [13, are not applicable in the present case. 

In order to understand this discontinuity, let us first consider a similar case in Cartesian coordinates as shown 
in Fig. m^a), where a TM wave Hi = ye^^"^+'^^-'^ is obliquely incident {k^ 7^ 0) from free space onto an uniaxial 
medium with permittivity tensor CxXX + Ctyy + ^tzz and permeability /it. The dispersion relation in this medium 
is k^^l {uP'tt^t) + k'^/ (ui'^e^^t) = 1. Thus, when is very small, kt^ becomes imaginary and the transmitted wave 
becomes evanescent. In the limit ~* 0, it can be found that the transmitted wave is so strongly evanescent that no 
fields exist in the region a; > 0. More interestingly, in the limit Cx ^ 0, the integration Etxdx has a finite value 
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FIG. 2: Reflection and transmission of (a) a TM wave in Cartesian coordinate system incident from free space onto an uniaxial 
medium whose e^c goes to zero; and (b) an outgoing TM wave in spherical coordinate system from the homogeneous and 
isotropic hole through a homogeneous uniaxial background medium whose radial permittivity goes to zero. 



(a) PEC (b) Special Uniaxial Medium 




FIG. 3: (Color online) Reflection of a right-handed circularly polarized wave from (a) a PEC and (b) a medium with normal 
permittivity and permeability go to zero simultaneously. The dotted arrows represent the direction of rotation. 

— 2irjQ COS Oe'^'^'^ /kz. In other words, Etx is compressed on the interface like a delta function. This special finite and 
nonzero value can be named as an electric surface voltage Ve- When a free charge q moves to the interface, this voltage 
will push it to the other side and transfer energy qVe to it. This voltage is not caused by conductive charges but an 
infinite polarization of the material on the interface, i.e. it corresponds to a distribution of polarized dipole moments 
on the interface. In addition, the tangential electric field at the left side of the interface is Ei^ + E^z = — St^q cosOe^'''^ 
while that at the right side is zero, meaning the tangential E field is discontinuous across the interface. However, 
since {Eiz + Erz)^z + Ve{z2) — Ve{zi) = 0, as shown in Fig. [Ha), Faraday's law still holds on this interface. Clearly, 
using this uniaxial material, which is the same with the inner boundary of the cloak, E^ becomes a delta function and 
forms the electric surface voltage which supports the discontinuity of the tangential E field. Meanwhile, the reflection 
coefficient becomes -1, meaning that this medium behaves like a PMC by means of controlling the medium's electric 
response. 

Similarly, the reflection coefficient for a TE wave is also -1 if goes to zero. Thus, the special uniaxial medium 
whose ex and Hx go to zero simultaneously behaves like a PMC for TM waves due to electric surface voltages, and a 
PEC for TE waves due to magnetic surface voltages. This leads to another interesting aspect of the reflection. First, 
in the sense that there is a complete reflection, the interface behaves like a mirror. Second, this special mirror not 
only keeps the polarization but also the phase information of the reflected waves. For example, for a right-handed 
circularly polarized wave incident onto this boundary, the reflected wave retains its handncss, but for a mere PEC 
or PMC boundary, the reflected wave becomes left-handed, as shown in Fig. [3| This property is similar to the SHS 
boundary used in radar and microwave engineering p^ . But for a SHS with its conducting vector flxed, if the incident 
plane changes, the phase of the reflected wave also changes. However, the phase of the reflected wave in Fig. Elb) 
is independent on the incident plane, meaning it only depends on the optical path the wave travels. So, this mirror 
behaves the same in any plane of incidence, and the information of a source including the polarization and phase is 
entirely retained in the reflected wave. 

From the above discussion in Cartesian coordinate system, we see that the surface voltages are introduced by the 
zero permittivity and permeability in the normal direction of the interface, which contribute to the discontinuity of the 
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FIG. 4: (Color online) Distributions of (a)the amplitude of Ve at the inner boundary and (b) Hy in the concealed region 
r < Ri due to an unit electric dipole oriented in z direction and located at {Ri/2, ■7r/4, tt) indicated by a small arrow in (b). 
Ri = 1.33Ai. 

tangential electromagnetic fields across the boundary. This is also true for a spherical interface in spherical coordinate 
system. For example, as shown in Fig. [2^b), a sphere with permittivity ei and permeability /ii is embedded in the 
homogeneous background medium with permittivity tensor e = e^ff + et09 + et(t>(f> and permeability fit- Similar to 
the case in Fig. [DJa), for TM waves, in the limit ~^ 0, no fields exist in the region r > Ri, but the electric surface 
voltage Ve is induced at the boundary. Since E™*Rid9 + Ve{6 + dO) — Ve{9) = 0, Faraday's law still holds across the 
boundary. The similar result can be obtained for component. Outgoing TE waves have similar derivation when 
fir goes to zero. Therefore the condition where the material at the inner boundary of a spherical cloak has radial 
permittivity and permeability of zero is sufficient for total reflection of all waves back by inducing surface voltages, 
no matter whether the outside medium satisfies the relation of constitutive parameters proposed in Ref. [l[ or not. 
Mathematical treatment in time domain in Ref. [l9| has also gotten the similar result of complete reflection. 

Based on the above discussion, the electric and magnetic surface voltages at the inner boundary as well as the 
field distribution inside the concealed region due to an electric dipole p located at an arbitrary position (r ,6 ^(j) ), 
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where r < Ri, can be derived. By expanding the wave from the dipole into spherical waves, the corresponding scalar 
potentials '^^m ^^'^ ^te f*-*^ incident waves can be obtained. Since it is known that the reflection coefficient for 
both TE and TM waves is ~(n{kiRi)/ipn{kiRi), the scalar potentials of reflected waves, '^^e ^^'^ ^tmj '^^^ be easily 
obtained. Consequently, the induced electric and magnetic surface voltages at the inner boundary of the spherical 
cloak can be calculated as follows 

Ve = r Erdr = — + ^Wl)\r^R^ , (8) 

Jr- ^^A^i^i cir 1 



Hrdr = — , 

R- t^Miei or 



Vh - / Hrdr = + ■^TE)\r=RT ■ (9) 



Figure 3] plots the amplitude of Ve at the inner boundary of a spherical cloak and the field Hy inside the concealed 
region in the x-z plane, due to an electric dipole pointing in z direction and located at (i?i/2, 7r/4, tt). Firstly, it 
is seen that surface voltages distributions are not uniform on the surface. But for an outside observer, the dipole is 
invisible since no wave propagates outside. Secondly, the field inside exists in the form of standing waves. Figure l^b) 
shows the field at the moment that the magnetic field reaches maximum. After a quarter of cycle, the magnetic field 
becomes zero while the electric field reaches maximum. Since E and H are always out of phase, the time-averaged 
Poynting power is zero anywhere, meaning no time-averaged power flowing inside. In other words, the energy radiated 
from the dipole at this moment will be returned to the dipole the next moment. Thus the total energy inside will not 
blow up. 

It can be calculated from Eqs.® and ^ that, in the presence of an electric dipole inside, when lo decreases to 
zero, Vh becomes zero while Ve survives. Similarly, if a static magnetic dipole is inside instead of a static electric 
dipole, Ve vanishes while Vh survives. The cloak for the static magnetic field can be realized artificially [13] ■ Since 
there is no magnetic charge in nature, this magnetic surface voltage induced by a static magnetic dipole must exist in 
the form of its equivalent electric surface current. The inner boundary in Ref. (20| is made of superconductor which 
makes this surface current realizable. 

Furthermore, the value of surface voltages can relate to another parameter directly. In derivation of the scalar 
potential in Ref. [ll|, the condition ^\E'tm = «we/x(^e [1^1 , where ipf. represents the auxiliary electric scalar potential, 
is applied. It is interesting to note that Ve = ^e\r=R-- Similarly, Vh = Vm\r=R- ■ Thus these auxiliary scalar 
potentials, (pe and (prm which were introduced originally as mathematical tools, have direct physical counterparts at 
the inner boundary of the cloak, i.e. surface voltages in this case. 

In conclusion, the exact electromagnetic field solution to the Maxwell equations for a spherical cloak with an active 
source inside the concealed region is established. It is shown that, an infinite polarization of the material at the inner 
boundary of the cloak will induce the electric and magnetic surface voltages, which prevent all waves from going out. 
These peculiar surface voltages are rare in nature, but they do not violate the Maxwell equations. The handness of 
the polarization and phase information of the waves reflected from the inner boundary of the cloak are unchanged. 
Finally, these surface voltages due to an electric dipole inside the concealed region are found to be exactly equal to 
the auxiliary scalar potentials at the inner boundary of the cloak which gain physical counterparts in this case. 

This work is supported by the ONR under Contract N00014-01-1-0713, the Chinese NSF under Grant 60531020, 
and China PSF under Grant 20060390331. 
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